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/.  Introduction 


We  have  used  rf  teehniques  with  plasma  probes  in  laboratory  experiments  to 
demonstrate  the  existenee  of  eollisionless  resistanee  in  the  sheath  of  a  spherieal  probe', 
shown  that  this  leads  to  a  method  of  finding  the  eleetron  sheath  density  profile^,  and 
proposed  a  method  of  measuring  eleetron  temperature.^  In  the  most  reeent  work  we 
demonstrated  a  method  of  determining  plasma  potential."'  In  the  present  work  we  show 
that  the  eleetron  distribution  funetion  ean  be  determined  from  a  1*^'  derivative  of  the 
inverse  ae  resistance  as  covered  below.  The  usual  convention  requires  a  derivative  of 
collected  electron  current  with  respect  to  applied  voltage  bias  in  the  neighborhood  of  the 
plasma  potential  since  it  is  only  there  that  electron  and  total  probe  collected  current  are 
nearly  the  same. 


11.  Determining  f(s) 


Druyvesteyn^  showed  that  the  electron  energy  distribution  function,  f(e),  for  a 


negatively  biased  probe  in  a  plasma  is  given  by, 

4  ms 


m= 


^  d^i  ^ 


V  P  JeV=-e 


(1) 


where  Ap  is  the  probe  area,  4  is  the  electron  current  collected  by  the  probe,  me  is  electron 
mass,  and  Vp  is  the  bias  voltage  with  respect  to  plasma  potential,  (pp,  and  e  is  the  energy. 

Total  current  collected  by  a  Langmuir  probe.  Ip,  will  include  both  electron  current 
in  addition  to  any  negative  ion  current,  along  with  positive  ion  current.  Because  of  this 
d  le/dVp  can  only  be  set  equal  to  d  Ip/dVp  for  small  biases,  Vp.  This  implies  that  f(s)  can 
be  estimated  for  low  energies  but  not  high  ones.  Common  procedure  relies  on  assuming 
a  form  for  the  ion  current,  which  can  be  subtracted  from  the  total  to  give  4.  As  this  form 
for  Ip  is  often  arbitrary  and  not  well  justified,  we  have  outlined  in  recent  work^  a  method 
of  obtaining  more  accurate  estimates. 


A  better  solution  would  involve  measuring  4  directly  without  the  necessity  of 
obtaining  it  from  the  total  collected  current  thus  eliminating  ion  current  contributions 
altogether.  As  discussed  above  when  applying  a  small-amplitude  rf  signal  to  a  probe,  the 
ions  contribute  little  to  the  total  current  if  the  applied  frequency,  co,  is  such  that  copi  «  a> 
where  cOpi  is  the  ion  plasma  frequency.  In  addition  if  co  <  C0pe{ro)  with  O0pe{ro)  the  electron 
plasma  frequency  at  the  probe  surface,  tq.,  we  avoid  electron  resonance  effects  discussed 
in  an  earlier  paper',  and  f(s)  can  be  expressed  as 


m  = 


dKl 

dV. 


\ 

JeV^ 


=-s 


(2) 


where  Rac  is  the  ac  resistance. 


2  2  -1 

The  justification  for  replacing  d  h/dVp  with  dRad  /dVp 
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is  covered  in  the  Appendix. 

HI.  Brief  Description  of  Experimental  Procedure 

We  refer  the  reader  interested  in  the  experimental  details  to  earlier  works'’"^  and 
only  provide  an  outline  of  that  same  description  here. 

The  experiments  we  deseribe  here  were  condueted  using  rwo  small  stainless  steel 
spheres  of  1.25  em  and  2.54  em  radius  whieh  were  eonneeted  to  an  HP8735D  Network 
Analyzer  through  50  Q  eoaxial  eable  whieh  provides  the  driving  signal.  This 
arrangement  ineluding  the  ehamber,  analyzer  and  the  eoupling  eireuitry  is  shown 
sehematieally  in  Figures  1  and  2.  The  spheres  were  mounted  on  a  1/4  ineh  diameter 
eeramic  and  steel  support  whieh  is  eonneeted  to  1/4  ineh  diameter  semi-rigid  eopper  50 
Ohm  coaxial  cable. 


For  all  of  the  experiments,  the  determination  of  plasma  impedanee  depends  upon 
the  network  analyzer  measurement  of  the  eomplex  refleetion  eoeffieient,  F(0).  From  this 
measurement  the  analyzer  returns  as  separate  outputs  Re  Zafoo)  and  Im  Zafoo)  where, 


l  +  r((u) 

i-r(®) 


(3) 


and  Zo  (=50  O)  is  the  internal  impedanee  of  the  analyzer, 
refleeted-to-total  power  is  given  by, 

Po 


We  also  note  that  the  ratio  of 


(4) 


where  Po  =  Pr  +  Pt  with  Pr  and  Pr  the  refleeted  and  transmitted  powers,  respeetively. 
(The  quantity  1-  |  F  |  ^  is  the  normalized  transmitted  power  and  this  output  is  also 
available). 

IV.  Results  for  f(s)  and  Comparison  to  Langmuir  Probe  Data 


Using  Eq.  (2)  we  plot  as  examples  f(e)  for  two  probes  of  radius  1.25  em  (Figure 
4),  and  2.54  cm  (Figure  5).  For  each  of  the  plots  we  have  determined  the  average  energy, 
Eav,  and  the  eleetron  density,  Ue,  from  the  usual  integrals  over  energy  of  ef(e)  and  f(e), 
respeetively.  The  eleetron  temperature  is  estimated  as  ~  2/3  Eav  The  plots  eaeh  show 
good  correspondenee  between  the  eonventional  Langmuir  probe  results  and  those  of  the 
method  presented  here. 
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V.  Summary 

We  have  presented  a  method  of  determining  the  eleetron  distribution  funetion 
using  teehniques  developed  over  the  past  three  years  in  studies  of  rf  impedanee  probes 
using  very  small  amplitude  signals  eompared  to  bias  levels  or  typical  plasma  potential 
magnitudes.  For  this  reason  our  technique  is  non-perturbative  to  the  existing  plasma- 
probe  interface  and  is  unaffected  by  the  presence  of  probe  contamination. 

VL  Appendix 


A.l  The  origin  ofEq.  (2) 


The  total  ac  impedance,  Zac,  of  a  circuit  whose  model  includes  a  parallel 
combination  of  Cac,Rac  to  represent  the  sheath,  in  series  with  a  bulk  plasma  inductance, 
Lp,  is  given  by. 


ZA^)  =  Z^^^{G))  +  ia)L^ 


1 


icoC„„  +  i?. 


-  +  i(oL^ 


(A.l) 


where  Zacjco)  is  used  to  denote  sheath  impedance,  Cac  is  sheath  ac  capacitance  discussed 
below  and. 


L 


P 


Anr^e^n^ir) 


(A.2) 


with  rs  denoting  sheath  radius.  (In  earlier  work  ,  we  included  a  sheath  inductance  but  this 
term,  which  considerably  complicates  the  solution,  was  found  to  be  negligible  for  the 
frequencies  we  consider.)  After  some  algebra  it  can  be  shown  that. 


R 


ac 


Re(Z.,) 


(A.3) 


and. 


= 


-Im(Z^J 


coicoL  -  Im(Z^^ ))  +  ®  Re(Z^^ ) 


Eliminating  Lp  from  Eqs.  (A.3)  and  (A.4)  we  are  able  to  write 


(A.4) 


where. 


Re(Z.,)  = 


R^ 


(A.5) 
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(A.6) 


Kc  =  ( 


die 

dV„ 


1 

/Id 

2 

4nsQ 

v'i)  y 

\ 

Te 

and  we  require  that  copi  «  co  <  copeiro)  with  copeiro)  the  eleetron  plasma  frequeney  at  the 
probe  surfaee,  rg.  Vp  is  taken  as  the  de  bias  with  respeet  to  plasma  potential,  Cac  is  sheath 
eapaeitanee,  Xd  is  Debye  length,  and  me,Te  are  eleetron  mass  and  temperature, 
respeetively.  Eqs.  (A. 3)  and  (A. 4)  arise  from  the  eireuit  representation  of  the  probe 
plasma  interaetion  whereas  Eq.  (A.6)  is  based  on  physieal  theory  whieh  ineludes  an 
assumption  of  Boltzmann  eleetrons  and  ignores  ion  eurrent  eontributions.^  Elnlike  the 
ease  where  a  very  low  frequeney  is  applied  {oo  <  oopi,  cOpe),  the  lower  bound  to  the 
frequeney  range  above  avoids  most  ion  eontributions  to  the  total  ae  eurrent.  For  the 
upper  bound,  there  will  be  no  eontribution  from  eollisionless  resistanee  (CR)  (or 
resonanee  effeets)  eovered  in  the  original  work  in  this  series.^  Eq.  (A.6)  and  the 
reasoning  whieh  produees  it  are  the  basis  of  the  transition  from  the  Druystevn  equation 
seen  in  Eq.  (4)  to  the  form  in  Eq  (5). 

If  Eq.  (A. 5)  is  solved  for  Rac,  we  obtain. 


\-^\-4co^Rc(ZJ^Cl 

Ico^ClRciZJ 


(A.7) 


And  so,  to  the  extent  that  {coRQ(Zac)Cacf  «  1  in  Eq  (A. 5)  ,  Rac  ~  RefZ^c) 
Therefore  Rac  has  a  minimum  at  cpp.  This  will  generally  hold  true  for  small  biases  near  cpp 
and  the  frequeneies  and  sheath  sizes  applieable  in  our  experiment,  i.e.,  frequeneies  no 
larger  than  approximately  O.lco^e,  and  probe  radii  no  larger  than  2.54  em. 


We  then  have  aiVp  =  (pp. 


and  finally. 


dV„ 


■  =  -R: 


v<y 


=  0 


dR^iZJ 


dV„ 


-0 


(A.8) 


(A.9) 


Beeause  of  this  we  ean  plot  the  network  analyzer  output  of  Re(Zac)  versus  applied 
bias  for  frequeneies  in  the  range  speeified.  Eaeh  of  these  eurves  will  then  show  a 
minimum  at  Fp  =  (pp.  An  example  of  these  plots  from  previous  work"^  for  the  two  probes  is 
shown  in  Fig.  (3).  Due  to  the  dependenee  of f(s)  on  dRac^/dVp  seen  in  Eq  (2)  we  are  able 
to  eonstruet/fe)  in  the  vieinity  of  cpp  for  small  negative  bias  voltages.  This  is  the  basis  of 
Figs  (4)  and  (5) 
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A,2  Determining  Cac 

Eq  (A. 4)  which  relies  almost  entirely  on  experimental  data  ean  be  used  to 
determine  ac  capacitance,  Cac,  from  Re(Zac),  Im(Zac),  and  Lp  as  found  from  Eq  (A. 2) 
where  ne(r)  is  determined  from  solution  of  the  Poisson  equation.  Note  that  the  lower 
integration  boundary  is  the  sheath  edge,  We  make  a  distinetion  here  between  Cac  and 
dc  capacitance,  Cdc-  This  ean  be  seen  when  realizing  that  we  are  applying  a  small  rf 
signal  to  a  biased  probe.  The  total  eharge  on  the  probe  may  then  be  expressed, 

0(0 = a. + ay"'  (A.io) 


where  Qac«  Qdc  and  therefore  to  lowest  order. 


Q  =V 

^ac  at 


dQ. 


dc 


dV^ 


dc 


with,  C  ^  = 


dc 


a 

E.. 


=  Anr^  a, 


2„  dE^{r^) 


0  ^0 


dV^ 


dc 


(A.  11) 


Vdc  is  the  potential  with  respeet  to  the  bulk  plasma  and  Vac  «  E/c  is  the  ae  voltage.  The 
expression  for  Cac  in  Eq.  (A.  11)  arises  from  the  faet  that,  for  a  spherieal  probe  in  a  bulk 
plasma. 


=  4m-ie,Edr,)  and,  (A.12) 

Eq  (A.12)  results  from  observing  that  the  eleetrie  field  vanishes  far  from  the  probe  and 
therefore  the  probe  and  the  bulk  plasma  eontain  equal  amounts  of  eharge.  Earlier  work^ 
ean  provide  Er(ro),  the  eleetrie  field  at  the  surfaee,  as  a  funetion  of  the  applied  bias.  From 
this  model,  it  is  possible  to  numerieally  fit  Er(ro)  as  a  funetion  of  Vdc  and  to  determine 
dEr(ro)/dVdc  and  henee  Cac-  For  a  typieal  ease  we  find  that  Cac  is  similar  to  Cdc  with  Cac 
being  generally  larger.  However,  the  magnitudes  of  those  differenees  do  not  substantially 
alter  the  results  here,  i.e.,  the  effeet  on  the  determination  of  plasma  potential  and /(h)  is 
insignifieant  for  the  frequeney  ranges,  probe  sizes,  and  magnitudes  of  Re(Zac)  near 
plasma  potential.  This  ean  be  seen  in  the  denominator  of  Eq.  (A. 5)  where  the  term 
((oRacCacf  «  1  is  typieal  for  these  parameter  regimes.  Also,  eomparison  to  Langmuir 
probe  measurements  is  largely  eonsistent  with  the  plasma  parameters  (Ne,  Te)  ealeulated 
from  the  determination  of f(e)  using  the  method  outlined. 

It  is  also  useful  to  note  that  Cdc  can  be  be  estimated  in  two  ways  from  vaeuum 
values.  Using  either  as  an  estimate  of  Cac  does  not  greatly  affeet  the  results  of  the 
determinations  off(e).  However  sinee  both  are  estimates  of  Cdc  they  are  not  justified  on 
physieal  grounds  as  a  eorreet  representation  of  Cac-  One  method  uses  the  vaeuum 
eapaeitanee  of  an  isolated  sphere,  Co,  whereas  another  uses  the  vaeuum  eapaeitanee  of 
two  eoneentrie  spheres,  the  outer  defined  by  the  sheath  edge.  The  vaeuum  estimates  are 
also  useful  as  benehmarks  and  for  eomparison  to  the  results  of  Eq.  (A. 7),  i.e.,  using  the 
numerieal  model  for  a  sphere  of  radius  IOAdc  for  example,  yields  a  typieal  ratio  2.5  < 
Cdc/Co  <3.7  when  only  applying  a  de  bias.  In  this  ease,  the  lower  bound  eorresponds  to 
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plasma  potential,  Vp,  and  the  upper,  near  1  V  removed  from  Vp.  In  the  seeond  ease  is 
estimated  from  a  solution  of  the  Poisson  equation.  This  solution  proeeeds  by  uniquely 
determining  bulk  eleetron  density,  no  ,  from  the  resonanee  in  Re(Zac)  and  the  zero 
erossing  of  Im(Zac).  With  this  information  we  are  able  to  estimate  with  only  an 
assumption  as  to  the  approximate  eleetron  temperature  Te.  This  is  possible  beeause  the 
sheath  width  is  largely  independent  of  Te  for  this  model. ^  The  Poisson  equation  is  solved 
by  speeifying  these  boundary  eonditions  at  the  sheath  edge  and  integrating  inward  to  the 
final  bias  whieh  produees  a  given  and  therefore  the  two-sphere  vaeuum  eapaeitanee. 
This  method  produees  an  estimate  which  is  always  larger  than  that  of  an  isolated  sphere 
but  is  increasingly  unphysical  when  approaching  the  limit  of  a  vanishing  sheath.  In  an 
intermediate  range  of  bias  voltages  (i.e.,  for  sufficiently  large  sheaths)  this  method  is 
closer  to  the  results  of  the  correct  calculation  of  Cdc  and  varies  as  a  function  of  bias. 

An  advantage  to  having  the  analytic  expression  for  Cac  of  Eq.  (A.  11)  is  that  it  is 
possible  compare  to  experimental-based  results  obtained  through  Eq  (A.4)  to  a  more 
realistic  theoretical  basis  than  the  vacuum  estimations  described  in  the  last  paragraph. 
This  can  be  done  for  example  by  using  the  Poisson  equation-based  calculation  described 
above  to  find  Er(ro)  and  hence  by  differentiation,  Cac-  This  is  planned  in  future  work 
which  will  address  the  issue  in  more  detail. 


VIL  Figure  Captions 

Figure  (1)  -  A  photograph  of  the  present  NRE  Space  Physics  Simulation  Chamber 
showing  magnetic  field  coils  and  two  separate  experimental  areas  separated  by  a  large 
gate  valve  providing  either  experimental  coupling  between  the  two  or  isolation  for 
separate  experimentation.  The  results  presented  in  this  work  were  performed  in  the  larger 
section  surrounded  by  five  magnetic  field  coils  toward  the  rear  of  the  photograph. 

Figure  (2)  -  A  schematic  representation  of  the  Network  Analyzer  and  coupling  circuitry 
necessary  for  the  swept  frequency  analyses  of  the  spherical  impedance  probes.  The 
Analyzer  returns  a  representation  of  the  signal  reflected  from  the  plasma  and  provides 
plasma  complex  impedance.  The  circuitry  shown  indicates  both  the  application  of  the 
small  non-perturbative  analyzer  signal  in  addition  to  the  dc  bias  applied  to  the  probes. 

Figure  (3)  -  A  plot  of  the  Re(Zac)  for  frequency  scans  varying  from  11  to  20  MHz 
returned  by  the  Network  Analyzer  vs  applied  bias  voltage  for  the  25.4  mm  radius  sphere. 
For  this  run  the  plasma  potential  determined  from  Eangmuir  probe  data,  ~  1.5  V, 
plasma  temperature,  Te  ~  0.4  eV  and  electron  density,  Ue  ~  10^  cm'^. 


Figure  (4)  -  A  plot  of  the  electron  distribution  function/(h)  vs  applied  probe  bias  for  the 
2.54  cm  probe  for  an  analyzer  sweep  frequency  of  3  MHz.  From  f(s)  we  find  by 

o  ■y 

integration  that  He  =  1.23  x  10  cm'  and  Te  =  0.38  eV .  This  is  to  be  compared  to  a 

o  ^ 

Eangmuir  probe  determination  of  He  -  1  x  10  cm'  and  Te~  0.4  eV 
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Figure  (5)  -  A  plot  of  the  electron  distribution  function  f(s)  vs  applied  probe  bias  for  the 
1.25  cm  probe  for  an  analyzer  sweep  frequency  of  4  MHz.  From  f(s)  we  find  by 
integration  that  He  =  3.4  x  10  cm'  and  Te  =  0.74  eV  .  This  is  to  be  compared  to  a 
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Langmuir  probe  determination  of  He  -  4  x  10  cm'  and  Te~  0.8  eV 
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